3.3 \(\int (e x)^m (a+b x^2) (A+B x^2) (c+d x^2) \, dx\)

Optimal. Leaf size=97 \[ \frac{(e x)^{m+3} (a A d+a B c+A b c)}{e^3 (m+3)}+\frac{(e x)^{m+5} (a B d+A b d+b B c)}{e^5 (m+5)}+\frac{a A c (e x)^{m+1}}{e (m+1)}+\frac{b B d (e x)^{m+7}}{e^7 (m+7)} \]

[Out]

(a*A*c*(e*x)^(1 + m))/(e*(1 + m)) + ((A*b*c + a*B*c + a*A*d)*(e*x)^(3 + m))/(e^3*(3 + m)) + ((b*B*c + A*b*d +
a*B*d)*(e*x)^(5 + m))/(e^5*(5 + m)) + (b*B*d*(e*x)^(7 + m))/(e^7*(7 + m))

________________________________________________________________________________________

Rubi [A]  time = 0.0605376, antiderivative size = 97, normalized size of antiderivative = 1., number of steps used = 2, number of rules used = 1, integrand size = 27, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.037, Rules used = {570} \[ \frac{(e x)^{m+3} (a A d+a B c+A b c)}{e^3 (m+3)}+\frac{(e x)^{m+5} (a B d+A b d+b B c)}{e^5 (m+5)}+\frac{a A c (e x)^{m+1}}{e (m+1)}+\frac{b B d (e x)^{m+7}}{e^7 (m+7)} \]

Antiderivative was successfully verified.

[In]

Int[(e*x)^m*(a + b*x^2)*(A + B*x^2)*(c + d*x^2),x]

[Out]

(a*A*c*(e*x)^(1 + m))/(e*(1 + m)) + ((A*b*c + a*B*c + a*A*d)*(e*x)^(3 + m))/(e^3*(3 + m)) + ((b*B*c + A*b*d +
a*B*d)*(e*x)^(5 + m))/(e^5*(5 + m)) + (b*B*d*(e*x)^(7 + m))/(e^7*(7 + m))

Rule 570

Int[((g_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_.)*((c_) + (d_.)*(x_)^(n_))^(q_.)*((e_) + (f_.)*(x_)^(n_))^
(r_.), x_Symbol] :> Int[ExpandIntegrand[(g*x)^m*(a + b*x^n)^p*(c + d*x^n)^q*(e + f*x^n)^r, x], x] /; FreeQ[{a,
 b, c, d, e, f, g, m, n}, x] && IGtQ[p, -2] && IGtQ[q, 0] && IGtQ[r, 0]

Rubi steps

\begin{align*} \int (e x)^m \left (a+b x^2\right ) \left (A+B x^2\right ) \left (c+d x^2\right ) \, dx &=\int \left (a A c (e x)^m+\frac{(A b c+a B c+a A d) (e x)^{2+m}}{e^2}+\frac{(b B c+A b d+a B d) (e x)^{4+m}}{e^4}+\frac{b B d (e x)^{6+m}}{e^6}\right ) \, dx\\ &=\frac{a A c (e x)^{1+m}}{e (1+m)}+\frac{(A b c+a B c+a A d) (e x)^{3+m}}{e^3 (3+m)}+\frac{(b B c+A b d+a B d) (e x)^{5+m}}{e^5 (5+m)}+\frac{b B d (e x)^{7+m}}{e^7 (7+m)}\\ \end{align*}

Mathematica [A]  time = 0.0600721, size = 73, normalized size = 0.75 \[ x (e x)^m \left (\frac{x^4 (a B d+A b d+b B c)}{m+5}+\frac{x^2 (a A d+a B c+A b c)}{m+3}+\frac{a A c}{m+1}+\frac{b B d x^6}{m+7}\right ) \]

Antiderivative was successfully verified.

[In]

Integrate[(e*x)^m*(a + b*x^2)*(A + B*x^2)*(c + d*x^2),x]

[Out]

x*(e*x)^m*((a*A*c)/(1 + m) + ((A*b*c + a*B*c + a*A*d)*x^2)/(3 + m) + ((b*B*c + A*b*d + a*B*d)*x^4)/(5 + m) + (
b*B*d*x^6)/(7 + m))

________________________________________________________________________________________

Maple [B]  time = 0.006, size = 321, normalized size = 3.3 \begin{align*}{\frac{ \left ( Bbd{m}^{3}{x}^{6}+9\,Bbd{m}^{2}{x}^{6}+Abd{m}^{3}{x}^{4}+Bad{m}^{3}{x}^{4}+Bbc{m}^{3}{x}^{4}+23\,Bbdm{x}^{6}+11\,Abd{m}^{2}{x}^{4}+11\,Bad{m}^{2}{x}^{4}+11\,Bbc{m}^{2}{x}^{4}+15\,Bbd{x}^{6}+Aad{m}^{3}{x}^{2}+Abc{m}^{3}{x}^{2}+31\,Abdm{x}^{4}+Bac{m}^{3}{x}^{2}+31\,Badm{x}^{4}+31\,Bbcm{x}^{4}+13\,Aad{m}^{2}{x}^{2}+13\,Abc{m}^{2}{x}^{2}+21\,Abd{x}^{4}+13\,Bac{m}^{2}{x}^{2}+21\,Bad{x}^{4}+21\,Bbc{x}^{4}+Aac{m}^{3}+47\,Aadm{x}^{2}+47\,Abcm{x}^{2}+47\,Bacm{x}^{2}+15\,Aac{m}^{2}+35\,Aad{x}^{2}+35\,Abc{x}^{2}+35\,Bac{x}^{2}+71\,Aacm+105\,Aac \right ) x \left ( ex \right ) ^{m}}{ \left ( 7+m \right ) \left ( 5+m \right ) \left ( 3+m \right ) \left ( 1+m \right ) }} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((e*x)^m*(b*x^2+a)*(B*x^2+A)*(d*x^2+c),x)

[Out]

x*(B*b*d*m^3*x^6+9*B*b*d*m^2*x^6+A*b*d*m^3*x^4+B*a*d*m^3*x^4+B*b*c*m^3*x^4+23*B*b*d*m*x^6+11*A*b*d*m^2*x^4+11*
B*a*d*m^2*x^4+11*B*b*c*m^2*x^4+15*B*b*d*x^6+A*a*d*m^3*x^2+A*b*c*m^3*x^2+31*A*b*d*m*x^4+B*a*c*m^3*x^2+31*B*a*d*
m*x^4+31*B*b*c*m*x^4+13*A*a*d*m^2*x^2+13*A*b*c*m^2*x^2+21*A*b*d*x^4+13*B*a*c*m^2*x^2+21*B*a*d*x^4+21*B*b*c*x^4
+A*a*c*m^3+47*A*a*d*m*x^2+47*A*b*c*m*x^2+47*B*a*c*m*x^2+15*A*a*c*m^2+35*A*a*d*x^2+35*A*b*c*x^2+35*B*a*c*x^2+71
*A*a*c*m+105*A*a*c)*(e*x)^m/(7+m)/(5+m)/(3+m)/(1+m)

________________________________________________________________________________________

Maxima [F(-2)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Exception raised: ValueError} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((e*x)^m*(b*x^2+a)*(B*x^2+A)*(d*x^2+c),x, algorithm="maxima")

[Out]

Exception raised: ValueError

________________________________________________________________________________________

Fricas [B]  time = 1.50248, size = 566, normalized size = 5.84 \begin{align*} \frac{{\left ({\left (B b d m^{3} + 9 \, B b d m^{2} + 23 \, B b d m + 15 \, B b d\right )} x^{7} +{\left ({\left (B b c +{\left (B a + A b\right )} d\right )} m^{3} + 21 \, B b c + 11 \,{\left (B b c +{\left (B a + A b\right )} d\right )} m^{2} + 21 \,{\left (B a + A b\right )} d + 31 \,{\left (B b c +{\left (B a + A b\right )} d\right )} m\right )} x^{5} +{\left ({\left (A a d +{\left (B a + A b\right )} c\right )} m^{3} + 35 \, A a d + 13 \,{\left (A a d +{\left (B a + A b\right )} c\right )} m^{2} + 35 \,{\left (B a + A b\right )} c + 47 \,{\left (A a d +{\left (B a + A b\right )} c\right )} m\right )} x^{3} +{\left (A a c m^{3} + 15 \, A a c m^{2} + 71 \, A a c m + 105 \, A a c\right )} x\right )} \left (e x\right )^{m}}{m^{4} + 16 \, m^{3} + 86 \, m^{2} + 176 \, m + 105} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((e*x)^m*(b*x^2+a)*(B*x^2+A)*(d*x^2+c),x, algorithm="fricas")

[Out]

((B*b*d*m^3 + 9*B*b*d*m^2 + 23*B*b*d*m + 15*B*b*d)*x^7 + ((B*b*c + (B*a + A*b)*d)*m^3 + 21*B*b*c + 11*(B*b*c +
 (B*a + A*b)*d)*m^2 + 21*(B*a + A*b)*d + 31*(B*b*c + (B*a + A*b)*d)*m)*x^5 + ((A*a*d + (B*a + A*b)*c)*m^3 + 35
*A*a*d + 13*(A*a*d + (B*a + A*b)*c)*m^2 + 35*(B*a + A*b)*c + 47*(A*a*d + (B*a + A*b)*c)*m)*x^3 + (A*a*c*m^3 +
15*A*a*c*m^2 + 71*A*a*c*m + 105*A*a*c)*x)*(e*x)^m/(m^4 + 16*m^3 + 86*m^2 + 176*m + 105)

________________________________________________________________________________________

Sympy [A]  time = 1.93092, size = 1515, normalized size = 15.62 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((e*x)**m*(b*x**2+a)*(B*x**2+A)*(d*x**2+c),x)

[Out]

Piecewise(((-A*a*c/(6*x**6) - A*a*d/(4*x**4) - A*b*c/(4*x**4) - A*b*d/(2*x**2) - B*a*c/(4*x**4) - B*a*d/(2*x**
2) - B*b*c/(2*x**2) + B*b*d*log(x))/e**7, Eq(m, -7)), ((-A*a*c/(4*x**4) - A*a*d/(2*x**2) - A*b*c/(2*x**2) + A*
b*d*log(x) - B*a*c/(2*x**2) + B*a*d*log(x) + B*b*c*log(x) + B*b*d*x**2/2)/e**5, Eq(m, -5)), ((-A*a*c/(2*x**2)
+ A*a*d*log(x) + A*b*c*log(x) + A*b*d*x**2/2 + B*a*c*log(x) + B*a*d*x**2/2 + B*b*c*x**2/2 + B*b*d*x**4/4)/e**3
, Eq(m, -3)), ((A*a*c*log(x) + A*a*d*x**2/2 + A*b*c*x**2/2 + A*b*d*x**4/4 + B*a*c*x**2/2 + B*a*d*x**4/4 + B*b*
c*x**4/4 + B*b*d*x**6/6)/e, Eq(m, -1)), (A*a*c*e**m*m**3*x*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 15*
A*a*c*e**m*m**2*x*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 71*A*a*c*e**m*m*x*x**m/(m**4 + 16*m**3 + 86*
m**2 + 176*m + 105) + 105*A*a*c*e**m*x*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + A*a*d*e**m*m**3*x**3*x*
*m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 13*A*a*d*e**m*m**2*x**3*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m +
 105) + 47*A*a*d*e**m*m*x**3*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 35*A*a*d*e**m*x**3*x**m/(m**4 + 1
6*m**3 + 86*m**2 + 176*m + 105) + A*b*c*e**m*m**3*x**3*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 13*A*b*
c*e**m*m**2*x**3*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 47*A*b*c*e**m*m*x**3*x**m/(m**4 + 16*m**3 + 8
6*m**2 + 176*m + 105) + 35*A*b*c*e**m*x**3*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + A*b*d*e**m*m**3*x**
5*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 11*A*b*d*e**m*m**2*x**5*x**m/(m**4 + 16*m**3 + 86*m**2 + 176
*m + 105) + 31*A*b*d*e**m*m*x**5*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 21*A*b*d*e**m*x**5*x**m/(m**4
 + 16*m**3 + 86*m**2 + 176*m + 105) + B*a*c*e**m*m**3*x**3*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 13*
B*a*c*e**m*m**2*x**3*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 47*B*a*c*e**m*m*x**3*x**m/(m**4 + 16*m**3
 + 86*m**2 + 176*m + 105) + 35*B*a*c*e**m*x**3*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + B*a*d*e**m*m**3
*x**5*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 11*B*a*d*e**m*m**2*x**5*x**m/(m**4 + 16*m**3 + 86*m**2 +
 176*m + 105) + 31*B*a*d*e**m*m*x**5*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 21*B*a*d*e**m*x**5*x**m/(
m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + B*b*c*e**m*m**3*x**5*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) +
 11*B*b*c*e**m*m**2*x**5*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 31*B*b*c*e**m*m*x**5*x**m/(m**4 + 16*
m**3 + 86*m**2 + 176*m + 105) + 21*B*b*c*e**m*x**5*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + B*b*d*e**m*
m**3*x**7*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 9*B*b*d*e**m*m**2*x**7*x**m/(m**4 + 16*m**3 + 86*m**
2 + 176*m + 105) + 23*B*b*d*e**m*m*x**7*x**m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105) + 15*B*b*d*e**m*x**7*x**
m/(m**4 + 16*m**3 + 86*m**2 + 176*m + 105), True))

________________________________________________________________________________________

Giac [B]  time = 1.16127, size = 645, normalized size = 6.65 \begin{align*} \frac{B b d m^{3} x^{7} x^{m} e^{m} + 9 \, B b d m^{2} x^{7} x^{m} e^{m} + B b c m^{3} x^{5} x^{m} e^{m} + B a d m^{3} x^{5} x^{m} e^{m} + A b d m^{3} x^{5} x^{m} e^{m} + 23 \, B b d m x^{7} x^{m} e^{m} + 11 \, B b c m^{2} x^{5} x^{m} e^{m} + 11 \, B a d m^{2} x^{5} x^{m} e^{m} + 11 \, A b d m^{2} x^{5} x^{m} e^{m} + 15 \, B b d x^{7} x^{m} e^{m} + B a c m^{3} x^{3} x^{m} e^{m} + A b c m^{3} x^{3} x^{m} e^{m} + A a d m^{3} x^{3} x^{m} e^{m} + 31 \, B b c m x^{5} x^{m} e^{m} + 31 \, B a d m x^{5} x^{m} e^{m} + 31 \, A b d m x^{5} x^{m} e^{m} + 13 \, B a c m^{2} x^{3} x^{m} e^{m} + 13 \, A b c m^{2} x^{3} x^{m} e^{m} + 13 \, A a d m^{2} x^{3} x^{m} e^{m} + 21 \, B b c x^{5} x^{m} e^{m} + 21 \, B a d x^{5} x^{m} e^{m} + 21 \, A b d x^{5} x^{m} e^{m} + A a c m^{3} x x^{m} e^{m} + 47 \, B a c m x^{3} x^{m} e^{m} + 47 \, A b c m x^{3} x^{m} e^{m} + 47 \, A a d m x^{3} x^{m} e^{m} + 15 \, A a c m^{2} x x^{m} e^{m} + 35 \, B a c x^{3} x^{m} e^{m} + 35 \, A b c x^{3} x^{m} e^{m} + 35 \, A a d x^{3} x^{m} e^{m} + 71 \, A a c m x x^{m} e^{m} + 105 \, A a c x x^{m} e^{m}}{m^{4} + 16 \, m^{3} + 86 \, m^{2} + 176 \, m + 105} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((e*x)^m*(b*x^2+a)*(B*x^2+A)*(d*x^2+c),x, algorithm="giac")

[Out]

(B*b*d*m^3*x^7*x^m*e^m + 9*B*b*d*m^2*x^7*x^m*e^m + B*b*c*m^3*x^5*x^m*e^m + B*a*d*m^3*x^5*x^m*e^m + A*b*d*m^3*x
^5*x^m*e^m + 23*B*b*d*m*x^7*x^m*e^m + 11*B*b*c*m^2*x^5*x^m*e^m + 11*B*a*d*m^2*x^5*x^m*e^m + 11*A*b*d*m^2*x^5*x
^m*e^m + 15*B*b*d*x^7*x^m*e^m + B*a*c*m^3*x^3*x^m*e^m + A*b*c*m^3*x^3*x^m*e^m + A*a*d*m^3*x^3*x^m*e^m + 31*B*b
*c*m*x^5*x^m*e^m + 31*B*a*d*m*x^5*x^m*e^m + 31*A*b*d*m*x^5*x^m*e^m + 13*B*a*c*m^2*x^3*x^m*e^m + 13*A*b*c*m^2*x
^3*x^m*e^m + 13*A*a*d*m^2*x^3*x^m*e^m + 21*B*b*c*x^5*x^m*e^m + 21*B*a*d*x^5*x^m*e^m + 21*A*b*d*x^5*x^m*e^m + A
*a*c*m^3*x*x^m*e^m + 47*B*a*c*m*x^3*x^m*e^m + 47*A*b*c*m*x^3*x^m*e^m + 47*A*a*d*m*x^3*x^m*e^m + 15*A*a*c*m^2*x
*x^m*e^m + 35*B*a*c*x^3*x^m*e^m + 35*A*b*c*x^3*x^m*e^m + 35*A*a*d*x^3*x^m*e^m + 71*A*a*c*m*x*x^m*e^m + 105*A*a
*c*x*x^m*e^m)/(m^4 + 16*m^3 + 86*m^2 + 176*m + 105)